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ABSTRACT 

One of the most outstanding problems in the gravitational collapse scenario of early structure 
formation is the cooling of primordial gas to allow for small mass objects to form. As the neu- 
tral primordial gas is a poor radiator at temperatures T < 10"* K, molecular hydrogen is needed 
for further cooling down to temperatures T ~ 100 K. The formation of molecular hydrogen is 
catalyzed by the presence of free electrons, which could be provided by the ionization due to 
an early population of cosmic rays. In order to investigate this possibility we developed a code 
to study the effects of ionizing cosmic rays on the thermal and chemical evolution of primor- 
dial gas. We found that cosmic rays can provide enough free electrons needed for the forma- 
tion of molecular hydrogen, and therefore can increase the cooling ability of such primordial 
gas under following conditions: A dissociating photon flux with F < 10"'^ erg cm^^Hz"' s"'. 



initial temperature of the gas 
sources with ecu > lO"-''' 



erg cm s ' . 



10 K, total gas number densities n > 1 cm , Cosmic ray 
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1 INTRODUCTION 

In the commonly adopted theory for structure formation, primor- 
dial density perturbations grow through gravitational Jeans' insta- 
bility and eventually form halos in whose potential wells the ob- 
servable baryonic matter gathers. During contraction, the gas de- 
velops shocks and gets reheated to a temperature at which pressure 
support can prevent further collapse. Once the gas has virialized 
in the potential wells of pre-existing dark matter halos, additional 
cooling is then required to further collapse the gas and form lumi- 
nous objects. 

Line cooling, in which the radiative de-excitation of atoms 
leads to the emission of photons, is the most important cool- 
ing mechanism for a primordial gas (predominantly made of 
hydrogen and helium) at temperatures T > 10'' K. At lower 
t emperatures though, primordial gas is a very poor radiator 
Jshapiro & Kanglii987h and in the absence of elements with lower 
excit ation energies, such as metals, line cooling becomes ineffec- 
tive jOalli & Pallalll993) . In the standard cosmological hierarchi- 
cal scenario for structure formation the objects which form first 
are predicted to have masses corresponding to virial temperatures 
r,.,> < lO"* K (for a review see ICiardi & Ferrarall2005h . Thus these 
objects, once virialized, cannot cool and further collapse via atomic 
line cooling. 

The most efficient coolant for a primordial gas at such temper- 
atures is molecular hydrogen, because it has additional rotational 
and vibrational degrees of freedom and therefore is able to cool the 
gas down to temperatures T ~ 100 K. For this reason molecular 
hydrogen is believed to play an important role in the formation of 
stars or small-mass galaxies from a metal-free gas. 



Therefore, it is of great interest to get a firm picture of the pri- 
mordial chemistry and of all the effects that might lead to formation 
or destruction of molecular hydrogen in metal-free gas. At redshifts 
z < 110 molecular hydrogen is predominantly formed by the gas 
phase reactions 

H+H H2 + e" 

in which electrons act as a catalyst ( lDalgamo&Lepdll987h . Thus, 
the formation of molecular hydrogen requires free electrons to be 
abundant at relatively low temperatures (T < 10'* K). 

In a post recombination universe the gas is mainly neutral with 
a residual fraction of free electro ns which give rise t o the formation 
of traces of molecular hydrogen jOalli & Fallal 1998h . Nevertheless, 
this primordial abundance i s not enough to trigger structure forma- 
tion iHirata & Padmanabhanl2006l:lAbel et aI.ri997l:lTegmark et al.l 



1 19971 : 1 Haiman & Loeb|[l997h . Partial ionization of hydrogen and 
helium increases the abundance of free electrons and promotes the 
formation of molecular hydrogen. While UV photons can fully ion- 
ize such species and, in the absence of neutral hydrogen, prevent the 
formation of H2, x-rays only produce a partial ionization and pro- 
mote molecular hydrogen formation. In the past, several authors 
have investigated the influence of x-rays on the H2 chemistry (e.g. 
Glover & Brand (2003); Haiman et al. (2000)). 

Recently much work has been published on the effect on pri- 
mordial chemistry of high energetic cosmic rays aris ing from the 
decay of massive dark m a tter pa rti cles (e.g.,S hchekin ov & VasilievI 
|2004 'Ripamo nti et"^ |2006|) . IShchekinov & VasilievI j2004h 
pointed out that ultra high energetic cosmic rays (UHECR), aris- 
ing from decaying super heavy dark matter particles with masses 
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Mx > 10'^ GeV may have a strong effect on the chemical evolu- 
tion of the medium. Those high energy cosmic rays are believed 
to be converted into UV photons through electromagnetic cascades 
which would partially ionize hydrogen and eventually lead to the 
formation of molecular hydrogen. The authors showed that in the 
presence of UHECR the early luminous objects are an order of 
magnitude less massive and form earlier in comparison to the pre- 
dictions of the standard recombination history. 

Cosmic rays (CR) themselves can be an alternative source of 
partial ionization. Unlike photons, cosmic ray particles can ionize 
neutral atoms many times as they move through the medium. In 
addition, in each ionizing collision the cosmic ray particles im- 
part only a sma ll fraction of their overall kinet i c energy to the 
ejecte d electron jMannheim & Schlickeisejl 19941 : [Spitzer & ScottI 
Il969h . Therefore cosmic rays are expected to heat the medium 
much less than photons with comparable energies. There are sev- 
eral suggestions about possible cosmic ray sources in the early 
universe. Decaying primordial black holes, topological defects, the 
very early supernovae, super massive particles or structure forma- 
tion shocks are some o f these (Biermann & Sigl 2001 ; Stanev 2004; 
IPfrommer et^ll2006l) . 

In this paper we will investigate whether an early population 
of cosmic rays can trigger molecular hydrogen formation and con- 
tribute to the process of small scale structure formation. This will 
be done by coupling in a self-consistent way the spectral evolution 
of a cosmic ray population with the chemical evolution of a pri- 
mordial gas. In Section|2]and[3]we describe the chemical evolution 
of the gas and the evolution of the cosmic ray spectrum respec- 
tively. In Section |4] and [5] we describe how cosmic rays interact 
with the ambient matter and by which reactions they influence the 
primordial chemistry. The method and the code developed during 
this work will be described in Section|6] Finally, Section[8]contains 
a discussion on the results which will be summarized in Section|9] 



2 PRIMORDIAL CHEMISTRY 

As discussed in the Introduction the basic idea of this work is to 
study the effect of cosmic rays on a gas consisting of different pri- 
mordial chemical species. Among others, we want to investigate the 
relevance of cosmic rays for early structure formation, which takes 
place in an almost metal free universe between redshift z ~ 40 and 
z ~ 10. We therefore consider only the chemistry of primordial el- 
ements, hydrogen, H, helium. He and deuterium, D. These three 
species interact via a complicated network of chemical reactions, 
which changes the overall chemic al state of the gas. 

According to ISpitzerl ( 1 1978 ) the velocity distributions of the 
gas particles in the density regime considered in this work (1 to 
10"* cm"^) are always close to a Maxwellian equilibrium distri- 
bution. This allows to describe the thermal state of this multi- 
component gas with just one single kinetic temperature T for all 
particles. 

Although the velocity distributions of the particle species are 
close to that in thermal equilibrium, the level population of the dif- 
ferent species does not correspond to thermal equilibrium. In the 
density regime considered here though, the coronal limit can be ap- 
plied. This allows to treat all particle species, including molecular 
hydrogen, in their ground states as long as only moderat e UV fields 
with fluxes up to ~ 10 erg cm - s ' Hz ' are present jAbel et akl 
ll997l : IShulllll97 8'). 

Given the above approximations, a simplified approach can be 
used to describe numerically the chemical and thermal evolution 



of the gas. The chemical evolution of a gas is governed by inelastic 
collisions that change the number densities of interacting species by 
creating or destroying chemical compounds. The equation which 
governs the number density evolution is most generally given as: 

'-^=Y^Y.'^%,iT)n,n,.Y^pX,n,-n'-^. (1) 

j k I 

Here the subscripts i,j,k and / denote the different particle species, 
and the quantities Kj^iT) and /IJ are the rates for the chemical re- 
actions that create or destroy particles of type ; ( reference to the 
used rates is given in Appendix iBll. a' and f}' are the stoichiometric 
coefficients that account for the number of particles of type ; that 
are created or destroyed in each reaction. The last term in equa- 
tion l[T]l describes the change in particle number density due to the 
evolution of the volume V. 

Whereas the temperature evolution of the gas is described by: 

dT 2(r-A) 2d\\iV 

— = -^^ T , (2) 

dt 3nkB 3 dt 

where n = is the total number density and T and A are the 
heating and cooling functions respectively (Appendix |AJ. The first 
term on the right hand side of equation (|2j describes the net thermal 
heat input to the medium, while the second term describes the work 
done by compression or expansion of the gas. Note that here we as- 
sume that the gas behaves as an ideal mono-atomic gas (Spitz^ 
Il978l) . Effects of thermal conduction are neglected in equation 
which is a ppropriate for the density and temperature range consid- 
ered here Jspitzeilll978h . 



3 COSMIC RAY SPECTRAL EVOLUTION 

It is well known that radiation like cosmic rays can ionize neutral 
atoms and deeply influence the chemical evolution of a medium. 
Since p rotons are the dominan t cosmic ray species at least in our 
Galaxy (*P frommer et al.l2006h , here we will only consider the cos- 
mic ray protons. This is a fair assumption as nucleosynthesis pre- 
dicts X 93% of all baryon particles in the universe to be protons. 
Therefore the dominant species emitted by some early cosmic ray 
accelerator should also be protons. However, a small fraction of 
helium might be present in the early cosmic radiation. These a- 
particles can simply be treate d as four protons a nd hence be ab- 
sorbed in the proton spectrum tenfilin et al.ll2006h . 

Since the ionization cross section of hydrogen has a maximum 
at proton impact energies a; 25.0 keV, it is necessary to find an 
adequate representation of the cosmic rays at those low energies. In 
addition, their original spectrum will be modified due to the energy 
losses which cosmic rays experience in the frequent interactions 
with particles of the ambient medium. As no measurement of low 
energetic cosmic rays is available and the cosmic ray ionization 
rate depends on the spectrum, we need to model the low energy 
spectrum and its time evolution theoretically. 

In the following we will define some variables useful for our 
calculations. It is convenient to introduce the dimensionless mo- 
mentum: 



where pcR is the momentum of the cosmic ray proton and m,, is the 
proton mass. In addition, we can express the kinetic energy Tp as: 

T„(p) = {^\+p^-\)(m„c^) (4) 
and the dimensionless velocity /? as: 
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(5) 



where v is the particle velocity. 

With this notation the spectral evolution equation of the cos- 
mic ray proton number density n(p, t) can be written as follows: 



dn{p, t) dp(p, t)n(p, t) 



dt 



dp 



Q(p,t)- 



n{p, t) 
r(p, t) ' 



(6) 



where the dot denotes the derivative with respect to time, Q{p, t) is 
a source function for the cosmic ray protons and T(p, t) is a typical 
loss time for cosmic ray particles from the volume under consider- 
ation. 

Under the assumption that the injected spectrum obeys a 
power law with some spectral index s in momentum space, the 
source function can be written as: 



Q(p,t)- 



nipC- k(s, po) 



p-'0(p-po). 



(7) 



where ecR is the change of total cosmic ray energy density with re- 
spect to time, k(s, po) is a normalization depending on the spectral 
index s, &(x) is the Heaviside function and po is a lower momen- 
tum cutoff of the injected spectrum. The value of k(s,po) can be 
calculated as (EnBlin et al. 2006): 



2 3-i- 



(8) 



where Bj(a, b) denotes the incomplete Beta-function, and i- > 2 is 
assumed. 



4 INTERACTION BETWEEN COSMIC RAYS AND 
MATTER 

4.1 Coulomb losses 

Here we are interested in the problem of the energy loss of cosmic 
rays in the interstellar and intergalactic gas. The basis of all elec- 
tromagnetic interactions is the Coulomb scattering between electric 
charges, which allows for transfer of kinetic energy as the proton 
moves through a free electron gas. According to iCouldl ( il972h . the 
total energy loss of a proton by Coulomb losses in a plasma is given 
by: 



i dT„(p) \ ^ Ane^n, 



In 



2mgC^I3p\ p' 

flOJ 



pi 



(9) 



Here, ojpi = ^jAne-n^lm^ is the plasma frequency and is the 
number density of free electrons. The result is independent of the 
mass nip of the incident proton as long as the recoil of the proton 
can be neglected, i.e. ynigC <k nipC, with y = 1/ - jfi- Thus 
equation l|9) is valid for y <sc mplnig. In a partially ionized medium 
the total loss of kinetic energy can simply be calculated by adding 
the losses due to excitation or ionization to the Coulomb losses as 
calculated here CGould. 19721 . 



4.2 lonisation losses 

Ionization energy losses are important for all charged particles, 
and for particles other than electrons and positrons they domi- 
nate over radiation energy losses at all but the highest energies. 
These ionization losses are governed by the Bethe-Bloch equation 



jGroom & KleiiibOOO l). which for moderately relativistic charged 
particles assumes the form: 

idTp(p)\ 



dt 



V Zn-Lo 

4ne* v-i 
> Zn. 



In 



Im^c^p^ 
IzHy) 



(10) 



where iiz is the number density of atomic species with elec- 
tron number Z, Iz is t he ioniz a tion p otential, and Lq is called 
the stopping number te ieglerl [1999). The quantity h(y) = 
-sjl + lynielmp + (mjmpy is a correction factor that accounts for 
the maximum kineti c energy that can be im parted to a free electron 
in a single collision i Groom & Kleinl2000l) . 

The density correction factor 6z accounts for the screening ef- 
fect of a medium which becomes polarized as highly relativistic 
particles move through it (Martin & Shaw. .2003,) . In this case the 
electromagnetic field of the ions may not be at the assumed free- 
space value, but is reduced by the dielectric constant of the medium 
(Ziegler 199^ . The density correctio n factor 6z is usually negligi- 
ble for gases dOroom & Kleinll200(]|) . but is given here for com- 
pleteness: 



2v - Dz 

2y-Dz + az{(yuz-y)/lnl0f^ 




(yi.z < y) 

(yo.z <y < yi(i)i) 

(y < yoz) 



Here, y = ln(p), Dz = I - 21n(SWp,//z) and yo.z, Ji,z, Oz, 
kz are empirical c onstants which characterize the atomic species 
JSternheime r" 1952) and are given in table[T]together with the values 
for the ionization potentials Iz of molecular hydrogen and helium. 
All hydrogen measurements were done with molecular hydrogen, 
but not much difference is expected when these values are adopted 
for atomic hydrogen. 

As already pointed out earlier, the energy range of cosmic 
ray protons most important for the primordial chemistry is of or- 
der keV. The Bethe-Bloch equation in the form l llOt is not valid 
at these low energies, and a correction term must be added to the 
square brackets of equation lll Ot. The so called Bloch correct ion is 
important for slow particles ( IMannheim & Schlickeiseilll994t) . and 
for protons it can be written as: 



1 / i la 

W;, = 2('P(l)-Re>I'(l + - 



(12) 



where l'(x) is the digamma function, the logarithmic derivative of 
the gamma function, and a = 1/137 is the fine structure constant 
(Ziegler 1999). 

With this corrections the Bethe-Bloch equation can be used for 
particles with velocities greater than the characteristic velocity of 
the medium 's electrons, which is the orbital velocity in the K-shell 
of the atom dMannheim & SchhckeiseJl994l) . For atomic hydrogen 
the orbital velocity can be derived from the Bohr atomic model and 
yields: 

2 

j8o,*« = ^ « 0.0073 . (13) 
he 

For protons with velocities larger than Porbn the Bethe-Bloch 
equation is in excellent agreeme nt with the experimental values 
dMannheim & Schlickeisej|l994l) . 

At even lower particle velocities Ginzburg and Sy- 
rovatskii provide a useful io nization loss rate for slow ions 
dGinzburg & Svrovatskii|[l964}) . which for protons can be written 
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Element 


Z 


h [ eV] 


yo,z 




"z 


kz 


He 


2 


24.6 


2.0 In 10 


3.0 In 10 


0.98 


4.11 


H2 


1 


13.6 


1.76 In 10 


3.0 In 10 


0.34 


5.01 



Table 1. Atomic data of He and H2 a.s given bylStemheimej <1952h . The 
hydrogen measurements were done with molecular hydrogen, but no signif- 
icant changes for atomic hydrogen ai'e expected within the accuracy needed 
for our purpose tEnBlin et al.2006.) . 



dt J 



2.11 X lO^y ^WzZ — 



(14) 



For proton veloci ties lower than P ortit this agrees well with the ex- 
perimental data of lWhaliiT^ ( Il958h . Then the total ionization energy 
loss formula is given by: 



( dT,(p) ] 



2z [Lq + Lgioch] 
2.ilx 10-2^62 2z«zZt 



iP > Par bill cx 
OS </?<,,*,,) • 



Note that the derivation of the cosmic ray energy losses due 
to ionization neglects the details of the atomic shell structure. In 
order to obtain accurate energy losses, the particle electron colli- 
sions should be considered with detailed treatment of each target 
electron's orbital bonding, and as the particle slows down, a she ll 
correction term should be added to equation JlOb (Ziegle3ri999h . 
All the atoms considered in this work though, on average will al- 
ways be in the ground state and thus no detailed study of the atomic 
shell corrections is needed. The c orrection wo uld be at most 6% 
in the energy range 1-100 MeV ( Izieglejri999 ). Therefore, shell 
correction can be safely neglected, and equation il5\ describes ad- 
equately the ionization losses of cosmic ray particles. 



4.3 Hadronic losses 

The known cosmic ray spectrum exte nds over energies from a 
few hundred MeV to 3 x 10^° eV (Bier mann & Siglll200l[) . Such 
high relativistic energies allow cosmic ray protons to overcome the 
Coulomb barrier and interact hadronically with the nuclei of the 
ambient matter The fundamental theory underlying high energy 
collisions of hadrons is generally referred to as Quantum Chromo 
Dynamics (QCD) (Stanev 2004). As QCD predicts the creation of 
massive particles in collisions between hadrons if the incident par - 
ticle possesses enough kinetic energy jPeskin & Schroededl 1995b . 
these interactions, especially pion creation, are an important energy 
loss mechanism for cosmic ray protons at high energies. The light- 
est known mesons are pions, with masses m„±c^ = 139.57 MeV 
for charged pions and m„oc^ = 134.98 MeV for neutral pions 
( lMartin& Shawl l2003h . For positive and neutral pion production 
the lightest final states which conserve baryon number and electric 
charge are created by the processes 



P + P-^P-i-N4-;t+ 
and 

P-nP^P-^P-F/r" 



(16) 



(17) 



yielding a threshold energy Z?"^ x 1.22 Ge V, with a correspondin g 
momentum threshold p,i,mpC^ x 0.79 GeV jMartin & Shawll2003h . 

The neutron produced in reaction ( |16b will decay into a proton 
after a mean lifetime of 886 s, so we basically end up with pions 
and two protons. This assures that the cosmic ray proton number 
density is conserved in these reactions. In subsequent processes the 
pions will decay via the channels 



^ p- + Vfj/v^ e- + vjve + v^/v^, 
and 





7t" 27 



(18) 



(19) 



with mean lifetimes fo r charged pions of 2.6 x 10"* s and for neutral 
pions of 8.4 x 10"" s jMartin & Shawll2003h . However, the subse- 
quent decay products are not considered further in this work. 

The average kinetic energy loss of cosmic ray protons due to 

J ion p roduction has been calculated by iMannheim & Schlickeiseil 
and can be written in the form 



l dT,(p) \ 
\ dt j„ 



■ 0.65cn„„,i,iCrpi,Ti,&(p - p,,,) , 



(20) 



where n,u,ciei is the number density of the target nucleon in the 
medium and cr,,,, = 3 x 10"^' cm^ is the cross section for pion 
production in proton collisions. 



5 COUPLING BETWEEN COSMIC RAY PROTONS AND 
PRIMORDIAL GAS CHEMISTRY 

In this Section we describe the method used and the approximations 
adopted to couple cosmic ray protons to the chemical evolution of 
a primordial gas. 

5.1 Ionization rates 

As a cosmic ray proton travels through matter it leaves behind 
a trail of ions. In accelerator experiments it is possible to mea- 
sure the mean energy needed to create one electron-ion pair Wz 
alon g the track, whic h clearly depe nds on the ta rget material 
l fBakker&Segrelll95lh . According to iHeitlej l fl954l) this average 
energy loss is representative for the average primary ionization of 
the cosmic ray proton. 

Both the fractional number of cases in which a collision with 
an atom results in ionization rather than excitation to a discrete 
level and the average energy transferred to the ionized electron, 
are al most independent of the incident energy of the cosmic ray 
proton jMannheim & SchUckeise3ll994h . With these empirical ob- 
servations and equation il5t , the primary ionization rate for the 
chemical elements can be calculated as: 



^liPmin) 



(21) 



Wztiz 

where £i(pmin) is the change of ionization energy density per time 
given as: 

dT„(p)\ 



£,( 



dt 



n(p)dp . 



(22) 



Here p„„„ is defined by the minimum energy necessary to ionize 
the atomic species (in the case of atomic hydrogen it is 13.6 eV) 
and Wz is the mean energy expended per i on pair created by a pro- 
ton. The values of Wz were m easured by iBakker & Segrel ( Il95 ih 
and lWeiss & BemsteinI ( Il956h for several elements, and are given 
in table |2] for molecular hydrogen and helium. The hydrogen mea- 
surements were done with molecular hydrogen, but no significant 
changes for atomic hydrogen are expected within the accuracy re- 
quired here. 

When a cosmic ray proton ionizes an atom, energy is trans- 
ferred to the ejected electron as kinetic energy. Some of this energy 
will go into further ionization or excitation of atoms, followed by 
radiation, while some of the energy will go directly into heating 
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the medium dSpitzer&Scotll 19691) . Therefore, this first generation 
electrons give rise to additional ionization and more free electrons 
as long as their kinetic energy is above the ionization threshold en- 
ergy. 

According to lSpitzer & Scotll il969h all secondary and higher 
generation ionization can be accounted for multiplying the primary 
ionization rate M\\ by a factor ^ = 5/3. This holds up to an ion- 
ization fraction of / = n^/nu « 10"^, while for higher ionization 
fractions secondary electrons loose their energy to Coulomb inter- 
actions with t he free electron gas rather than to ionization of addi- 
tional atoms iSpitzer & Scott 1969). As cosmic rays can raise the 
ionization fraction considerably, especially for low density gases, 
one should be aware of the errors that are introduced by using a 
constant ^. However, at this stage of research we are mainly in- 
terested in finding an overall effect and therefore we treat ^ as a 
parameter. Comparison of calculations with ^ = 5/3 and f = 1 
showed that not much difference arises from the different treat ment 
of secondary electrons in the cases studied here jjasc 

With all the above simplifications the total ionization rate, in- 
cluding primary and higher generation ionizations, can be written 
as: 



ifor J- ^ prill! 



(23) 



5.2 Cosmic ray input on the chemical network 

The most important reactions by which cosmic ray protons influ- 
ence the primordial chemical network are the following: 

(i) PcR + H ^ PcR + H+ + e- 

(ii) PcR -H He ^ PcR + He+ + e" 

(iii) PcR + D -» PcR + D+ -I- e- 

(iv) PcR + H2 ^ Pgr + HJ + e" 

(v) PcR -H H- ^ PcR -H H + e- 

where the electron detachment reaction (v) turned out to be negli- 
gible ( |jaschell2006h . 

All these reactions conserve the cosmic ray proton number 
density, and therefore no exchange of particles between cosmic rays 
and medium needs to be considered. The ionization rates for the 
reactions can be calculated with equation ([23} and the data taken 
from tables [T] and |2] Unfortunately, no experimental data for deu- 
terium targets is available, and therefore we use the same values as 
for hydrogen, which should be a safe assumption for the accuracy 
required in this work. As 97% of cosmic ray-Hi impacts lead to the 
formation of H^ via reaction (iv) no other destruction mechanism 
for molecular hydrogen by cosmic r ay protons will be considere d 
here in agreement with the literature jHartquist & Williamslll99^ . 
However, inclusion of additional cosmic ray-H2 reactions lead to 
the formation of additional H2. Such reactions and also the influ- 
ence of cosmic rays on other molecules like HD will be discussed 
in detail in forthcoming publications. 

Whenever there is helium present in the medium, it will not 
only be ionized by the cosmic ray protons, but will also emit a 
21.23 eV phot on, as the remaining electron relaxes to a new shell 
configuration jSpitzer& Scotlll969h . In principle this photon is ca- 
pable of ionizing further atomic hydrogen. But since we consider a 
very rarefied medium and the flux of these photons is assumed to 
be low, they are expected to contribute insignificantly to the overall 
ionization rate. 



Element 


Z 


Wz [ eV] 


He 


2 


40.3 ± 0.8 


Hz 


1 


36.3 ± 0.7 



Table 2. Experimental values of Wz for m olecular hydrogen 
iBakker & Segre 1951) and helium iWeiss & Bernstein. 1956) . 



5.3 Heating rates 

Cosmic rays can be an efficient heat source especially for a low 
density gas, and therefore their heat input to the medium must be 
considered. When the cosmic ray proton ionizes an atom, it trans- 
fers a certain amount of kinetic energy to the electron, which is 
either used for further excitation and ionization of atoms, or is dis- 
tributed by elastic collisions to other species of the medium. In the 
latter case the overall kinetic temperature is raised. 

The effect of heating by cosmic ray ionization and secondary 
electrons is not well described in the literature. The papers that 
mention heating by cosmic ray ionization are not very deta i led an d 
show some discrepanc ies (e.g. iMannheim & Schlickeiseil ( Il994l) : 
ISpitzer&Scot^ ( ll969l) ). Here we will provide a reasonable upper 
limit to the heat input. 

Let be the ionization threshold energy of an atom of 
atomic charge Z (for hydrogen Ej^^^ = 13.6 eV). Then the mean 
kinetic energy imparted to the primary electron in a single col- 
lision can be calculated as 



E„„ = {Wz-El) 



(24) 



In subsequent collisions the primary electron will on average ion- 
ize ^ - 1 atoms and therefore loose the additional kinetic en- 
ergy (<f - 1)^^,- Thus every primary ionization by a cosmic ray 
particle finally yields a temperature gain equivalent to an energy 
Eiiea = {Wz - ^Ef^i^- Hence we can write the heating rate for cos- 
mic ray ionization as 



J^IOII 

i CR 



:{Wz-^;EQA^'nz^ 



A-c'rIIz 



(25) 



which, for <f = 5/3, means a heat input of =a 8.18 eV for 
every ionization of hydrogen in a neutral medium. This heat in- 
put is an overestimation as the electron could have also lost energy 
to excitation which was not accounted for here. Thus, in a neutral 
medium less than 1 /4 of the energy transferred from the cosmic ray 
proton is finally deposited as heat. Again, it should be noted that f 
decreases with increasing ionization fraction / and thus the heat 
input by cosmic rays will increase with /. As soon as a large frac- 
tion of the medium is ionized, the cosmic ray proton looses kinetic 
energy mainly due to Coulomb interactions with the free electrons 
of the medium, and all the Idnetic energy imparted to the free elec- 
trons is equivalent to the heating energy Ei,i,„. 

The heating rate due to Coulomb interactions can be calcu- 
lated by using equation (|9}: 



T^COIll 

i CR ■ 



r 



dT„(p) 
dt 



n(p) dp . 



(26) 



The fact that the energy lost by cosmic ray protons in Coulomb in- 
teractions is directly transferred to heat of the medium is the reason 
why cosmic rays are an efficient heat source. Thus, once the gas is 
ionized to sufficiently high degrees the future fate of the medium 
will be defined by Coulomb heating solely. 
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6 METHOD 

Here we will describe and discuss the numerical methods used to 
study the influence of early cosmic rays on primordial gas. In order 
to estimate the elFects of cosmic ray protons on the chemistry of 
primordial gas, it is necessary to follow the chemical and thermal 
evolution of such gas along with the cosmic ray spectral evolution. 
Therefore we simulate a homogeneous and isotropic medium and 
follow the evolution of 13 species (H, H+, H , H:!^, H2, He, He^, 
He**, HeH*, D, D*, HD and e ) self-consistently together with the 
spectral evolution of the cosmic ray protons. The code calculates 
the chemical evolution of primordial gas by taking into account 
27 collisional rates ( Appendix IbI ). 17 radiative rates including 
CMB interactions ( Appendix |B] ) and 4 additional rates arising 
from the coupling of cosmic ray protons to the chemical network 
as mentioned earlier. 



6.1 Code testing 

The accuracy of the numerical scheme for the cosmic ray spectral 
evolution is checked by comparing the numerical results for the 
cosmic ray spectral evolution with the analytic equilibrium solution 
of equation l|6) where we neglect the escape term. The deviation of 
the numerical results from the analytic ones is measured by 



6n(p) 



(P) - nnumericalip) 

^analytic (P) 



(27) 



Here we do not consider any chemical evolution and assume a 
gas consisting of atomic hydrogen with a number density ;jh = 
10^ cm"^^. The cosmic ray spectrum is integrated from the lower 
boundary p„„„ = 10"^ to the upper boundary p„,nj- = 10^ with a 
spectral index i- = 2.2 and a cosmic ray energy density injection 
rate ecR = 10"^* ergcm"^ sec"'. 

The asymptotic steady state solutions are then compared for 
two scenarios. In the first scenario we consider that the cosmic 
rays interact with a totally neutral medium, while in the second 
we choose a totally ionized gas. The results are plotted in figure[T] 
The numerical solutions are nearly identical to the analytic ones. 
At the highest momenta the boundary effects lead to some devia- 
tions, which are still within the accuracy needed for this work. In 
any case these will not affect the calculations considered here since 
those high momenta are far beyond the region of interest for cosmic 
ray ionization. Therefore the calculations done in this work will not 
sufi^er from accuracy losses due to numerical effects. 

Testing the chemical network is more difficult as no analytic 
solution exists. For this reason, we have chosen to check the be- 
havior of the code against results present in the literature. We run 
a simulation to follow the chemical evolution of the intergalactic 
medium from redshift lO'' to and compare the results with those 
oflOaUi & Palla ( 1998). 

For our calculations we adopt a ACDM cosmological model 
with Six = 0.73, Sl„, = 0.27 and Ho = 71 kmsec"' Mpc"'. As 
initial conditions we assume a completely ionized gas with mass 
fractions /h = 0.76, /hc = 0.24 and /d = 1 .8 x 10"^ /h, as predicted 
by nucleosynthesis ( Smith etal. 1993). The result is presented in 
figure [2] From a comparison between the results obtained in this 
work and those of Galli and Palla, it is clear that the main features 
of the chemical evolution are very well reproduced by the code, 
although some minor differences are present. These arise because 
of the slightly different sets of chemical rates used and because 
Galli and Palla included more chemical species. 



Finally as chemical reactions conserve the mass of the system 
we checked that the mass remains constant throughout the calcula- 
tions. 



6.2 The simulation 

The gas in the simulations consists of the primordial elements H, 
D and He wit h the mass fraction s /h = 0.76, /hc = 0.24 and fa = 
1.8 X 10-^ /h (S mith etal JI1 993'). 

The parameters that determine the cosmic ray spectral evolu- 
tion are set as follows. We set the boundaries for the integration 
region in momentum space at p„„„ = 10"^ and p„,i,x = lO'', and 
assume the momentum cutoff po in equation ^ to coincide with 
p„,j„. The spectral index s of the injected power-law spectrum is as- 
sumed to be = 2.2, as expected for Fermi-acceleration at shocks 
jMannheim & SchlickeiseJl994h . This gives k(s, p„„,) = 4.95. The 
upper boundary is much higher than the range of momentum 
expected to be important for the chemical network, and therefore 
boundary effects arising from the numerical integration of the cos- 
mic ray spectrum will not interfere with the calculations of the 
chemical network. In order to restrict the number of parameters, 
the escape of cosmic ray particles from the medium under consid- 
eration is assumed to be negligible and therefore their escape time 
is set to T = 00. This implies that either the gas occupies an infinite 
volume, or magnetic fields capable of confining the cosmic rays in 
the region of interest are present. As already discussed in Section 
IS.ll in the following we will usually adopt ^ = 5/3 and a mean heat 
input per ionization of Ei,„„ = 8. 1 8 eV for hydrogen and 4.0 eV for 
helium. 

As not much about primordial cosmic ray sources and the re- 
gions in which cosmic rays might be present in the early universe 
is known, the code was used to run a large number of simulations 
in order to explore a wide range of the scr x " parameter space. 

Tests with an additional UV photon background field, as it 
might arise from the very first stars, showed that in the presence 
of ionizing photons the influence of cosmic rays is usually sub- 
dominant. We therefore consider only regions that are shielded 
from ionizing photon radiation, but the treatment of Ht dissoci- 
ating photons with energies in the ran ge of the Lyma n and Werner 
transitions (1 1.2-13.6 eV) is included jThielensl2005h . As only this 
small range of photon energy is of interest, we may assume a con- 
stant photon flux F. In the following we will refer to these dis- 
sociating photons as soft ultraviolet (SUV) radiation. No radiative 
transfer of photons is considered in this work, and therefore the 
dissociating background flux must be understood as a mean flux 
of photons in the gas. In reality, the molecules in the outer regions 
of the gas cloud absorb some of the dissociating photons resulting 
in a decrease of the dissociating radiation toward the inner regions 
( lThielensll2003) . This self-shielding mechanism is clearly depen- 
dent on the gas density and distribution and therefore the mean flux 
of photons in the medium will in general depend on the total num- 
ber density n and the sources that produce it. 



7 RESULTS 

In this Section we will discuss the effects of CR protons on the 
chemistry of a primordial gas as derived by our simulations. 

It should be noted that cosmic rays are not able to destroy an 
initially available population of molecular hydrogen. As the elec- 
tron detachment from H" by cosmic rays is negligible the only pos- 
sible destruction mechanism for molecular hydrogen is H2 ioniza- 
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tion. Nevertheless, as soon as cosmic rays are strong enough to 
sufficiently ionize molecular hydrogen, they also provide enough 
free electrons b y ionizing th e atomic species, resulting in a net pro- 
duction of H2 (Jasche 2006). Therefore, cosmic rays do not inter- 
fere with any other chemical process that leads to the formation of 
molecular hydrogen. This suggests that the final abundance of H2 
can be written as: 

«H, «= "hj + A«H, (28) 

where is the initial number density of molecular hydrogen, and 
A«Ht is the number density produced by cosmic rays in a purely 
atomic gas of the same total number density. 

The simple dependence of the overall behavior of /ih, on its 
initial abundance allows to make some more general considera- 
tions. Let us assume to have a gas cloud with initial H2 abundance 
> in thermal and chemical equilibrium in the absence of cos- 
mic rays. If H2 line cooling is the main cooling mechanism, this 
implies that T = A^^{n^l^J, where F includes all heating sources 
with exception of cosmic rays. When cosmic rays interact with the 
gas cloud, its temperature change can be written as: 

jnkeAT = T + r„ - Ah, 

= T + TcR- [ Ah2 ("h, ) + ( AnH2 )] 

= TcR - An,(Ann,) , (29) 



if we assume that all additional heating sources remain constant 
and equal to the initial H2 -cooling. Here we have used the fact that 
Ahj in thermal and chemical equilibrium is proportional to the H2 
number density. This implies that the temperature change induced 
by the interaction of cosmic rays and a gas with any initial H2 abun- 
dance can be estimated simply from the balance between Tcr and 
the H2 cooling of an initially atomic gas. The advantage is that we 
can explore a wider physical condition range without knowing the 
exact initial abundance of H2 and the details of heating sources 
other than cosmic rays. 

Note that this separation of the chemical network into a purely 
atomic part, which is influenced by the cosmic rays, and an un- 
changed part with initial molecular hydrogen, must not be under- 
stood as a rule. This approximation is only valid under the assump- 
tions presented here, and if the number density of atomic hydrogen 
does not change much. 

Therefore, we will start our calculations from a neutral purely 
atomic gas of primordial composition, and all molecular hydrogen 
H2 will be understood as the additionally produced molecular hy- 
drogen Awht • 

7.1 H2 and temperature 

The behavior of the chemical network under the radiation by cos- 
mic ray protons is not well known. What we expect from the cou- 
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Figure 2. Chemical evolution of a primordial gas described in this work (left panels) and in lGalli&Palljh998h (right panels). 



pling between the primordial chemical network and cosmic rays is 
an enhancement in the number density of H2, as cosmic rays ion- 
ize the atomic species and provide additional free electrons. This 
would eventually result in a net cooling of the gas by H2 line emis- 
sion if: 

(i) H2 is produced on sufficiently short timescales and H2 line 
cooling always outweighs cosmic ray heating; 

(ii) the gas is ionized only moderately to provide sufficient 
atomic hydrogen allowing for high H2 formation rates. 

The line cooling function Ah, of molecular hydrogen depends on 
the gas temperature T. As in a low temperature gas the colliding 



particles on average do not possess enough kinetic energy to ex- 
cite the rotational modes of molecular hydrogen, cooling by H2 
line emissions becomes more and more inefficient at lower temper- 
atures. Therefore we expect a minimum temperature at which no 
further cooling by molecular hydrogen is possible and cooling is 
balanced by the heating of cosmic rays. 

As Ah2 depends on wh, > sufficient amounts of molecular hy- 
drogen must be produced on short timescales to radiate away the 
heat input to the medium. For the formation of H2, free elec- 
trons and atomic hydrogen are needed. As cosmic rays ionize the 
medium they produce additional free electrons, but at the same time 
reduce the number of neutral hydrogen atoms. Therefore we expect 
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Figure 3. The left panel shows the final temperature and the right panel displays the final H2 number density as a function of ecR for two different SUV 
background fluxes. Different curves refer to different total number densities: n = 1 (solid line), 10 (dashed line), 100 (dotted line) and 1000 cm"^ (dash-dotted 
line). The initial temperature is Tq = 300 K. 



a lower H2 formation rate once a considerable fraction of the atomic 
hydrogen is ionized. 

To assess all these effects we run simulations with gas of four 
different total number densities (1-1000 cm"') and cosmic ray en- 
ergy injection rates ecR ranging from 10"'* to 10"^* erg cm"' s"', 
with different SUV background fluxes F. The evolution of the 
chemical network is followed over a time of 0.47 Gyr during which 
£cR is assumed to be constant. After this integration time the gas has 
reached nearly thermal and chemical equilibrium, and not much 
change is observed beyond it. 

The final values for the temperature T and the number density 
of molecular hydrogen jjht are plotted in figure [3] and figure |4] as a 
function of ecR for simulations with initial temperature Tq = 300 K 
and 3000 K respectively. 

The overall observed effect is counterintuitive, as the gas cools 
with higher energy injection by cosmic rays. As can be seen in the 
plots, the temperature behavior (left panels) is highly dependent on 
the given environmental conditions, i.e. the total number density n, 
the initial gas temperature Tq, the cosmic ray source strength €cr 
and the SUV background flux. As expected from the previous dis- 
cussions the final temperature is no monotonic function of the cos- 
mic ray energy density injection rate ecr- Beyond a certain value of 
£cR, which also depends on the total gas number density n, cosmic 
rays start heating the gas again instead of inducing further cooling. 



In addition we observe that gas with lower density is more likely 
to be heated by cosmic rays than gas of higher density, as here pro- 
duction of molecular hydrogen is easier and H2-line cooling is more 
efficient. 

The data also show a considerable impact of the SUV back- 
ground on the thermal evolution of the medium. A strong dis- 
sociating background suppresses the formation of molecular hy- 
drogen and thus H2 line cooling counterbalances the cosmic 
ray heating less efficiently. With a dissociating flux of F = 
10"'* ergcm"^ s"' Hz"' the production of molecular hydrogen is 
suppressed by a factor of 10"' compared to the fiducial calculations 
with F = 10""' ergcm"^ s"' Hz"', resulting in effectively stronger 
heating by cosmic rays. 

The choice of the initial temperature is critical for the evo- 
lution of the gas, as can be seen by a comparison between figure |3] 
and figureO which immediately shows that the effect of cosmic ray 
induced cooling is much more efficient at high initial temperatures. 
This is due to several effects, the most important being the increase 
of the H" formation rate with temperature and the fact that H2-line 
cooling is more efficient at higher temperatures. Beside these ef- 
fects the slight decrease of the recombination rates of the species 
H, D and He with temperature helps to increase the population of 
free electrons needed to catalyze the formation of H2. 
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Figure 4. Same as figure|3]but for initial temperature Tq = 3000 K. 



It is interesting to note that withi increasing initial temper- 
ature To cosmic rays of lower intensity become more and more 
important. While, for example, for a gas with Tg = 300 K 
and n = 10 cm"^ the first observable cooling effect begins at 
ecR ~ 10"^° erg cm"' s"', the cooling effect starts already at ecr ~ 
10^^^ erg cm"' s"' for the same gas simulated with To = 3000 K. 
This suggests that in the latter case cooling induced by cosmic rays 
might be efficient also at number densities < 1 cm"'. 



A closer comparison of figure [3] and figure |4] reveals another 
interesting feature. At the lowest energy injection rates the produc- 
tion of molecular hydrogen in gas To = 3000 K is increased by 
nearly one order of magnitude compared to the calculations with 
gas starting from To = 300 K. In contrast, at the highest simulated 
values of €cr the results for temperature T and Mh, ^6 identical for 
both initial temperatures. In particular the results of calculations 
with these two initial temperatures are identical for values of ecR 
higher than the energy injection rate at the minima in the tempera- 
ture curves. Therefore beyond this point the chemical and thermal 
evolution of the gas does not depend on the initial conditions any 
more and all gas starting from whatever initial condition is con- 
verted to the same final state. Thus the minimum in the temperature 
curve defines a characteristic point for each gas. 



7.2 Time evolution of parameter space 

So far we have discussed the behavior of primordial gas under the 
influence of cosmic rays comparing only the final results of several 
calculations. It is instructive to study also the temporal evolution of 
the gas. Therefore, we display the evolution of the ecr x «-plane 
in time as a sequence of images. To have a satisfactory resolution 
of the parameter space we have run calculations for 40 different 
values of n, in the range 1-10* cm"', and of ^cr in the range 10"'*- 
lO"^"* erg cm"' s"' . Each image of the €cr x ;i-plane then consists of 
1600 points, where each point represents a different simulation. 

In figure |5] we present some snapshots of a simulation with 
To = 2000 K and F = 10"^' ergcm"^ s"' Hz"'. It can be seen that 
molecular hydrogen forms from the very begitming of the simula- 
tion, at times when no change in temperature can be observed yet. 
Initially we observe the highest effect in correspondence with high 
values of the energy density injection rate £cr and total gas number 
density n, while the effect for low values of ecr and n takes time to 
build up. 



8 DISCUSSION 

The calculations presented in this work give us a better understand- 
ing of the behavior of a primordial gas exposed to radiation of cos- 
mic rays. In particular the parameter space maps, presented in fig- 
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Figure 5. Time evolution of the ecR x "-plane for calculations with a neutral atomic gas at initial temperature To = 2000 K, embedded in a dissociating 
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ureO can be used to identify regions of interest in the early uni- 
verse. Here we want to discuss these results in more detail. 

The most interesting conclusion that can be drawn from our 
calculations is that cosmic rays do not necessarily heat a gas, quite 
the contrary they might contribute to cooling. However, this trend 
is non-monotonic, as once the cosmic ray flux becomes larger than 
a critical value that depends on the gas initial conditions, heating 
dominates. This indicates that, to estimate the correct behavior of 
the thermal evolution of a primordial gas under the influence of 
cosmic rays, H2 chemistry should be included self-consistently. 

We studied the response of primordial gas to cosmic ray radia- 
tion in a variety of different environmental settings. In particular we 
were interested in the influence of the total gas number density n, 
the initial temperature To, the SUV background and the cosmic ray 
intensity represented by ecR. We observed that the effect of cosmic 
rays depends strongly on the environment, especially on the ini- 
tial temperature and the strength of the dissociating photon back- 
ground. The results presented here therefore enable us to clearly 
define the environmental conditions for which this effect might be 
of interest. The following conditions are necessary for the cooling 
induced by cosmic rays to be efficient (although some cooling can 
be observed also for lower temperatures and densities): 

(i) A dissociating photon flux with F < 10""* erg cm"^ Hz"' s"' . 

(ii) Initial temperature of the gas ~ 10' K. 

(iii) Total number densities of the gas n > 1 cm"'. 

(iv) Cosmic ray sources with ecR > 10"'' erg cm"' s"'. 

We can therefore think of a variety of possible astrophysical scenar- 
ios in which cosmic rays might efficiently influence the chemistry 
of a primordial gas. For example, the first two conditions are in 
excellent agreement with the environment that can be found at the 
epoch of first structure formation, when objects with T,,,, < 10'* K 
collapse. As the cosmic microwave background is already much 
too weak to ionize, the only possible source for ionizing photon ra- 
diation would be stars that have not formed in sufficient amounts at 
these times. During the process of structure formation shock heated 
gas can contribute to the ultraviolet background by thermal emis- 
sions. However, these fluxes are expected to be small, especially 
for small structures, and in case of the SUV b ackground do not ex - 
ceed values of F ~ 10"^' ergcm"^Hz"' s"' jMiniati et al]|2004h . 
In addition, X-rays can contribute to ionization. These high ener- 
getic photons arise from free-free interactions of the electrons in 
the shock heated gas. As soft X-rays do not penetrate deeply into 
clouds and the flux of hard X-rays is small, the ionization is ex- 
pected to be dominated by cosmic rays in the regimes considered in 
this work (Thielens 2005). It should be noted that, to have a proper 
estimate of the relative importance of different ionization sources, a 
self-consistent calculation should be done. A more thorough study 
on effects of CR on early structure formation is deferred to future 
work. 

The main goal of the calculations presented above was to get a 
better understanding of the cosmic ray chemistry and effects which 
affect the H2 formation. Therefore in all calculations no volume 
evolution was considered. This introduces some systematic errors 
when we want to apply the results to real cosmological situations. 
As the gas expands or contracts when it is heated or cooled, the 
total number density of the gas particles changes as a response to 
temperature change. Therefore, when the gas starts to cool, its num- 
ber density increases leading to a higher formation rate of H2. This 
effect might increase the overall H2 formation rate for a gas which 
experiences cooling induced by cosmic rays. To properly assess 
the impact of the volume evolution on our calculations, a detailed 



treatment of the gas collapse and the effects of cosmic ray pressure 
should be included. This is beyond the scope of this study and we 
defer a more thorough investigation to future work. 

The conclusions drawn in this paper rely on the presence of 
CR protons in the early universe, but our knowledge about primor- 
dial CR is very limited. As already mentioned in the Introduction 
some possible sources of primordial CR are decaying primordial 
black holes, structure formation shocks or the very first Super- 
novae. There is also the possibility that weakly interact ing massive 
partic les (WIMPS) decay into protons and anti-protons jPiehl et aU 
I1995I) . This would be an interesting source for cosmic ray protons 
as it is homogeneously distributed in space. Unfortunately the pro- 
ton and antiproton fluxes arising from the decaying WIMPS are 
expected to be too low to have any interesting impact on the gas 
cooling (Diehl et al. 1995). 

A more promising cosmic ray source in the absence of stars 
may be structure formation shocks. Cosmological shock waves 
form abundantly in the course of structure formation, both due to 
infalling pristine cosmic gas, which accretes onto filaments, sheets 
and halos, as well as due to superson ic flows associated with merg- 
ing structures ( |Pfrommeretalj2006l) . These shocks are able to dis- 
sipate gravitational energy associated with hierarchical clustering 
into thermal energy of the gas contained in dark matter halos. Be- 
side this, shocks are able to accelerate ions of the high energy tail 
of the Maxwellian velocity distribution of a therm al medium by 
diffusive shock acceleration jPfrommer e7^l2006l) . This acceler- 
ation process produces a cosmic ray population with a power-law 
distribution of the particle momenta. For this reason structure for- 
mation shocks seem to be plausible cosmic ray sources. As in ad- 
dition we do not expect stron g UV fluxes during early structure 
formation dMiniati et alj2004h our calculations could be applied to 
such a configuration. Numerical studies suggest that the mean en- 
ergy density injection rate of cosmic rays from structure formation 
shocks atz = 10 is €cr ~ 10"" erg cm"' s"', where this value is av- 
erage d over a co-moving volume of « (143 Mpc)' jPfrommer et al.l 
l2006h . Thus, the local value of ecr might be much higher. 

Once the very first stars appear, much more powerful cos- 
mic ray sources become available. These early massive stars (M ~ 
IOOMq ) form from me tal-free gas through molecular hydrogen 
cooling jAbel et al.ll2 002t) and are lik ely to end either as supernova 
or a black hole ( Wis e & Abelll2005l) . Supemovae are believ ed to 
be very efficient particle accelerator teiermann&Sigll|200 1) and 
therefore could easily increase the abundance of an early cosmic 
ray population. This could help to make the effect observed in our 
calculations more efficient and trigger the formation of second gen- 
eration stars. Unfortunately primordial stars produce a lot of UV ra- 
diation which might severely hamper the production of molecular 
hydrogen ( Wise & A bel 2005). Also, as already pointed out earlier, 
in the presence of ionizing photons the effect of cosmic ray ioniza- 
tion is usually sub-dominant, and if the medium is partly ionized 
cosmic rays will rather heat by Coulomb heating than induce H2- 
line cooling. One should also note that stars produce metals which 
are far more efficient coolants than molecular hydrogen. Therefore, 
in order to observe the effects described in this work, cosmic rays 
have to extend further into space than the pollution by metals. For 
these reasons we do not expect any significant effect of cosmic ray 
induced cooling in the vicinity of the sites of the very first stars and 
supernovae. 

If, however, cosmic rays would be able to travel further away 
than the ionization front and the metal polluted area built up by 
the progenitor, they might contribute to the formation of molecular 
hydrogen in the ambient neutral medium. Anyway, such a scenario 
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would require detailed knowledge about the propagation properties 
of cosmic rays and magnetic fields in the early universe, which we 
presently do not have. 



9 CONCLUSION 

In this work we studied the influence of an early cosmic ray popula- 
tion on the chemistry of primordial gas. We developed a computer 
code to couple self-consistently the chemistry of a gas of primordial 
composition with the spectral evolution of cosmic rays permeating 
this medium. The code was then used to study the response of the 
gas in a variety of environmental settings. The main results of these 
calculations can be summarized as follows: 

(i) Cosmic rays do not necessarily heat, quite the contrary they 
can contribute efficiently to the cooling of a medium. 

(ii) Cosmic rays influence the chemical state of a gas apprecia- 
bly. In particular they catalyze the formation of molecular hydro- 
gen. 

(iii) The impact of cosmic rays on the gas chemistry depends 
highly on the environmental conditions. 

However, the detailed processes of cosmic ray production and 
propagation in the early universe are yet unclear and further inves- 
tigation will be needed to completely understand a potential rele- 
vance of cosmic rays for early structure formation. 
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Cooling rate Reference 



Reaction 



Reference 



A/,„.(H) 

A,,-,„(He) 

A,,„,(He+) 

A,„„(H) 

A,on(He) 

A,™(He+) 

A,o„(He) 

A,,,(H+) 

A,,,(He+) 

A*(He+) 

A„,(He++) 

^H-diss 

formation 



diss 
'^Ht line 



Cen (1992a) and Black (1981) 
C_en (1992a) and Black (1981) 
]Cen (1992a) and Black (1981) 
Cen (1992a) and Shapiro & Kang 
Cen (1992a) and Shapiro & Kang 
_Cen ( 1 992a) and Shapiro & Kang 
Cen ( 1 992a) and Shapiro & Kang 
Black (1981) and Spitzer (1978^ 
Black (1981) and Spitzer (19re) 
Black (1981) and Spitzer (1978) 
Black (1981) and Spitzer (1978) 
Black (1981) 
Peebles (1971) 
Shap iro & Kang (1987) 
A^bel et al. ( 1997) 
Shapiro & Kang (1987) 
Galli & Palla (1998) 



(1987) 
(198T) 
(198T) 
(1987) 



Table Al. Cooling rates and the corresponding references. 



APPENDIX A: HEATING AND COOLING RATES 

The cooling rates used in this work are given in table |AT] together 
with the corresponding references. 

The kinetic energy of particles (i.e. electrons) produced by 
photo-ionizat ion and photo-di s sociat ion will be the heat sources. 
According to lShapiro & Kanj ( Il987l) the heating rate due to these 
reactions is given by: 



T = n, r (/IV -K,,) F(v)o-'(y)^, 



(Al) 



where «, and vj^ are the density and threshold energy of species 
i, respectively, F(y) and (t'{v) are the radiation flux and the cross 
section for the reaction, respectively. Note that no radiative transfer 
is considered in equation [Al] 



APPENDIX B: CHEMICAL RATES 

The chemical reactions implemented in our code are given in table 
IB 1 [ together with the references to the chemical rates. 

This paper has been typeset from a TgX/ KTgX file prepared by the 
author. 



H + e- ^ 
H+ + e- - 
He + e" - 
He+ + e- 
He+ + e- 
He++ + e" 
H + e- 
H + H- - 
H + H+ - 
HJ + H - 
H2 + H+ 
H2 + e" - 
H2 + H - 
H- + e- - 
H- + H - 
H- + H+ 
H- + H+ 
HJ + e- - 
HJ + H- 
D+ + e- - 
D + H+ - 
D+ + H - 
D+ + H2 
HD + H+ 
He + H+ - 
HeH+ + 1 
HeH+ + e 
H + y ^ 
He + y - 
He+ +7 
H- + y - 

Hj+y- 
Hj+y- 
H2+y- 
H2+y- 
H + ycMB - 
He + ycMB ■ 
He+ + ycMB 
+ ycMB 
H2 + ycMB - 
Hj+ycMB -»H + H+ 
H+ + ycMB ^ 2H+ + e- 
D + ycMB ^ D+ + e- 
HeH++ycMB ^ He + H+ 



H+ + 2e" Janev & Langer (1987) 

^H + y Abel et al. (1997) 

» He+ + 2e" Janev & Langer (1987) 

^ He + y Cen (1992b) and Aldrovandi & Peguignot (1973) 

He++ + 2e- AMDIS Database (1989) see Abel et al. (1997) 

— > He^ + y scaUng (e.g. see^sterbrock ^989)) 

> H- + y Abel et al. (1997) 
H2 + e- Abel et al. (1997), Shapi ro & Kang (1987) 
HJ + y Shapiro & Kang (1987) 

->H2+H+ Karp asetal. (1979) 

HJ + H Abel et al, (1997) 

2H + e- Donahue & Shull (1991) 

^ 3H Dove et al. (1986) 

H + 2e" Janev & Langer (1987) 

2H + e- Abel et a l. (1997) 

^ 2H Dalgamo & Lepp (1987) 

-*Hj+e- Abeletal. (1997) 

2H Abel et al. (1997) 

-> H2 + H Dalga mo & Lepp (1987) 

D + y Gain & Palla (1998^ 

^ D+ + H GalU & Palla ( 1 998) 

-> D + H+ Galli & Palla ( 1 998) 

HD + H+ Smith et al. (1982) 

^H2 + D+ Smith etal. (1982) 

— > HeH^ + y Roberge & Dalgarrra 

1 ^ He + H J Karpas et al. ( 1 979™ 

r ^ He + H Yousif & Mitchell (1989) 

H+ + e" Osterbrock (1974) 

» He+ + e" Osterbrock (1974) 

He++ + e" Osterbrock ( 1974) 

+ H + e- de Jong (1972) 

> Hj + e" Oneil & Rei nhardt (1978) 
^ H + H+ Osterbrock (1974) 
* 2H+ + e" Shapiro & Kang (1987) 
>H*^H + H Abeletal. (1997) 
>H + H Abeletal. (1997) 

H+ + e- GalU & Palla ( 1 998) 

> He+ + e- Gain & Palla ( 1 998) 

^He+++e- jjani& Palla (1998) 

♦ H + e- Galh & Palla ( 1 998) 

>Hj+e- GalU & Palla (1998J 
GalU & Palla (1998) 
GalU & Palla (1998J 
GalU & Palla (1998) 
GalU & Palla (1998) 



1982) 



Table Bl. Reactions and references for the different chemical rates. 
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